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Abstract 

The O'Connell process is a softened version (a geometric lifting with a parameter 
a > 0) of the noncolliding Brownian motion such that neighboring particles can change 
the order of positions in one dimension within the characteristic length a. This process 
is not determinantal. Under a special entrance law, however, Borodin and Corwin 
gave a Fredholm determinant expression for the expectation of an observable, which is a 
softening of an indicator of a particle position. We rewrite their integral kernel to a form 
similar to the correlation kernels of determinantal processes and show, if the number of 
particles is N, the rank of the matrix of the Fredholm determinant is N. Then we give 
a representation for the quantity by using an iV-particle system of complex Brownian 
motions (CBMs). The complex function, which gives the determinantal expression to 
the weight of CBM paths, is not entire, but in the combinatorial limit a — > it becomes 
an entire function providing conformal martingales and the CBM representation for 
the noncolliding Brownian motion is recovered. 

Keywords The O'Connell process • Noncolliding Brownian motion • Geometric lifting 
• Combinatorial limit • Fredholm determinants • Quantum Toda lattice • Whittaker 
functions • Macdonald processes • Complex Brownian motions 

1 Introduction 

1.1 Background 

A determinantal point process is a random ensemble of points in a space such that all 
correlation functions are given by determinants, whose matrix entries are values of a single 
continuous function called the correlation kernel [28, 26]. It can be generalized to space- 
time systems and if all spatio-temporal correlation functions are given by determinants, the 
process is also said to be determinantal [9, 15]. The noncolliding Brownian motion with 
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a finite number of particles iV is determinantal for all deterministic initial configurations 
£(•) = J2f=i firj(-)- (For each AcR, £(A) = j A ^(dx) gives the number of particles in A in 
the configuration £. Especially, for r G M, £({r}) denotes the number of particles located at 
the point r.) In particular, if the initial positions of particles {rj}^ =l are all distinct (i.e., 
for r£l, £({ r }) = 1 if t — rj, 1 < j < N, and £({?"}) = otherwise), the spatio-temporal 
correlation kernels is explicitly given as, for (x,x') G M 2 , (t,t') G [0, oo) 2 [16] 

N 

K%(t,x;t',x')=Y^ / dyp(t,x\r j )p(t , ,y\0)$l i (x' + iy) - l {t>t , )P (t - t' ,x\x') (1.1) 
j=1 Jr 

with 

r:€(W)=l, 

where i = \/— 1, p(t, y|x) denotes the transition probability density of the one-dimensional 
standard Brownian motion (BM) 

e -(x-y) 2 /2t 

p(t,y\x) = 7 =^l (t>0 ) + 5(x-y)l (t=O ), (x,y)GM 2 , f > 0, (1.3) 

v 27rc 

and is the indicator function of a condition ou; l( w ) = 1 if u> is satisfied and l( w ) = 
otherwise. The results are extended to the infinite-particle systems, in which the function 
(1.2) is regarded as the Weierstrass canonical product representation of an entire function 
[16]- 

O'Connell introduced an iV-component diffusion process, N > 2, which can be re- 
garded as a stochastic version of a quantum open Toda- lattice [22]. The Hamiltonian of 
the GL(iV, R)-quantum Toda lattice is given by 

H a N = --A+^V N (x/a), x = (x u x 2 ,...,x N ) eR N (1.4) 

with the Laplacian A = J2f=i d 2 /dx 2 and the potential 



JV-1 

V N (x) = J2 e ~ (Xj+1 ~ Xj) - ( L5 ) 

The Weyl chamber of type A N _i is given by W N = {x = (x 1: x 2 , ■ ■ . , x N ) G R. N : Xi < 

x 2 < ■ ■ ■ < x N }. For v = (ui, u 2 ,..., v N ) G W N , the eigenfunction problem H%ip^\x/a) = 

\(u)ip^\x / a) for the eigenvalue \(v) = — \v\ 2 /2 is uniquely solved under the condition 

that e- v - x / a ^\x/a) is bounded and lim e~ v - x ' a ^\x/a) = J[ T(u k -u 3 ), 

x^oo,xew N i<j<k<N 

where x — > oo, x G Wat means Xj + i — Xj — > oo, 1 < j < N — 1, and T denotes the Gamma 
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function. The eigenfunction ip v (■) is called the class-one Whittaker function [2, 23]. The 
infinitesimal generator of the O'Connell process is given by [22] 

= ^A + Vlog^Or/^-V, (1.6) 

where V = (d/dx\, . . . , O/Oxn). This multivariate diffusion process is an extension of a one- 
dimensional diffusion studied by Matsumoto and Yor [20, 21]. (The Matsumoto-Yor process 
describes time-evolution of the relative coordinate of the N = 2 case.) 

We can show that the O'Connell process is realized as the mutually killing BMs con- 
ditioned that all particles survive forever, if the particle position of the j-th BM is iden- 
tified with the j-th component of the O'Connell process, 1 < j < N [12, 13, 14]. Let 
Bj(t), 1 < j < iV be independent one-dimensional standard BMs started at Bj(0) = Xj G R, 
and for v = (u u u 2 , . . . , v N ) G R N , B^ 3 (t) = Bj(t) + Ujt, 1 < j < N be drifted BMs. We con- 
sider an iV-particle system of BMs with drift vector v, B v (t) = (B^(t), . . . , Btf (t)) , t > 0, 
such that the probability P w (t|{S I/ (s)} < s <t) that all N particles survive up to time t con- 
ditioned on a path {S I/ (s)} < s <t decays following the equation 

^-P a N (t\{B»(s)} 0<s<t ) = -\v N (B»(t)/a)PUt\{B»(s)} 0<s<t ), t > 0. (1.7) 
at a z 

It is a system of mutually killing BMs, in which the Toda- lattice potential (1.5) determines 
the decay rate of the survival probability depending on a configuration B u (t) [13]. With the 
initial condition B u (0) = xe W N , the survival probability Pft(t) = E x [P^(t\{B u (s)} < s <t)] 
is obtained by averaging over all paths of BMs started at x, and we can show that [14, 23, 13] 

lim P£(t) = c a 1 (N,iy)e- iy - x/a ^\x/a), if v G W N , u ^ 0, 

Pn(!) ~ c a 2 {N)r N ^ N - l ^ A ^\x/a) as t^ oo, if v = 0, (1.8) 

where c1(N,u) and c 2 {N) are independent of x and t. Then, conditionally on surviving of 
all N particles, the equivalence of this vicious BM, which has a killing term given by the 
Toda-lattice potential, with the O'Connell process is proved. We note that the parameter 
a > in the killing rate (1.7) with (1.5) indicates the characteristic range of interaction to kill 
neighboring particles as well as the characteristic length in which neighboring particles can 
exchange their order in R. It implies that if we take the limit a — > 0, the O'Connell process 
is reduced to the noncolliding BM. (The original vicious Brownian motion is a system of 
BMs such that if pair of particles collide they are annihilated immediately. The noncolliding 
BM is the vicious BM conditioned never to collide with each other, and thus all particles 
survive forever.) 

In the present paper, the limit a — > is called the combinatorial limit and an inverse 
of this procedure is said to be a geometric lifting in the sense of [4]. (See also [3].) Since 
determinantal functions associated with noncolliding diffusion processes (e.g., the Karlin- 
McGregor determinants, the Vandermonde determinants, the Schur functions) are replaced 
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by functional of the class-one Whittaker functions in the geometric lifting [2, 22, 23], the 
O'Connell process is not a determinantal process. 

Recently Borodin and Corwin [6] introduced the family of probability measures on se- 
quences of partitions, which are written in terms of the Macdonald symmetric functions 
and specified by the Macdonald parameters q,t e [0, 1) [19] as well as two Macdonald non- 
negative specializations [6]. This family of discrete measures is not determinantal. They 
showed, however, that if we consider a sub-family of processes with t — called the q- 
Whittaker measures, and if we observe a special class of quantities, which are eigenvalues 
of Macdonald's difference operators and called the Macdonald process observables [6], then 
some determinants appear, though still the processes are not determinantal. Taking a q — > 1 
limit of the q- Whittaker measures leads to a family of Whittaker measures which are now 
continuous and supported in R N . The determinants survive this limit transitions. The 
interesting and important fact is that the Whittaker measures are also realized as probabil- 
ity distributions of the O'Connell process started according to a special entrance law. Let 
X a {t) = (Xf(t),X$(t), . . .,Xfr(t)),t > be the O'Connell process with N particles. This 
special entrance law is realized by setting X a (Q) = —Mp with (2.14) and taking the limit 
M — > oo as explained in Section 2.3 below (see also [22, 10]). The drift vector of N particles 
of the O'Connell process, v = (vi, z/ 2 , . . . , vn) £ R N plays a role as the parameters specifying 
the Whittaker measures of Borodin and Corwin. We will write the expectation with respect 
to the O'Connell process with N particles and drift vector v started according to this special 
entrance law from "— oop" as E"'°[-]. 

For x e R, a > 0, set 

Q a (x) = exp(-e- x/a ). (1.9) 
l( x >o), that is, (1.9) is a softening of an indicator function l( x >o). 



Note that lim© a (a;) 

For a configuration £(•) = J2f=i ^j(-)> denotes a simple positively oriented contour on 

C containing the points {rj}^ =l located on R. Let 5 = sup{\vj\ : 1 < j < N} and choose 

< 5 < 1 so that 8 < 5/2. Borodin and Corwin [6] proved that E">°[e o (Xf (*) - h)}, heR 
is given by a Fredholm determinant of a kernel K eh / a for the contour integrals on C(— v), 

"(•) = E?=iM0; 



v,a 



e a (x»(t) - h) 



Det 



S(V — V') + K £ h/a (v, v') 



;i.io) 



where 



/ioo+8 r N 
— r(- s )r(i + s )n 
■ioo+<5 Z7r< f _.. 



r(v + v t ) u s e tvs / a2+ts2 i 2a2 



L r(s + V + V() v + s — v 



— , u>0. (1.11) 



Here the Fredholm determinant is defined by the sum of infinite series of multiple contour- 
integrals 



Det 



oo L p 

'5(v-v') + K u (v,v')] =J2jjUf 



L=0 j=l 



'C(-u) 



^-^JKM], (1.12) 
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where the term for L = is assumed to be 1. Note that (1.11) depends on u,a and t; 
K u (-, •) = K u (-, •; v,a,t). 

The Fredholm determinant formula (1.10) discovered by Borodin and Corwin [6] is sur- 
prising, since the O'Connell process is not determinantal as mentioned above. We would like 
to understand the origin of such a determinantal structure surviving in the geometric lifting 
from the noncolliding BM to the O'Connell process. 

It is well-known in quantum mechanics that the wave function of iV-particle system of 
free fermions is expressed by an N x N determinant called the Slater determinant. Then, a 
determinantal process is also called a fermion point process [26]. One should be careful, how- 
ever, that the notion of fermion is not enough to formulate determinantal processes, since in 
the context of stochastic processes repulsive interactions between paths in a spatio-temporal 
plane should be described. In a previous paper [17] , as an extension of notion of free fermions, 
we gave the complex Brownian motion (CBM) representation for the noncolliding BM. Let 
Zj(t),t > 0, 1 < j < N be a set of independent CBMs such that the real and imaginary 
parts, denoted by Vj(t) = $lZj(t),Wj(t) = ZsZj, are independent one-dimensional standard 
BMs. Since given by (1.2) is entire, (Z,-(f)) is a conformal map of a CBM, and hence 

it is a time change of a CBM. In other words, ^ (Zj(t)), 1 < j < N provide a set of indepen- 
dent complex local martingales, which are called conformal local martingales in Section V.2 
of [24]. Therefore a determinant of N x N matrix, det [$WZ jt (-))] ) is a martingale for 

l<j,k<N ? 

the system of independent CBMs. We proved that the noncolliding BM can be represented 
by the system of independent CBMs weighted by this determinantal martingale [17]. In 
comparison of the CBM representation for the noncolliding BM with a free fermion system, 
free-ness of particles is ensured by independence of CBMs, Zj(-), 1 < j < N and fermionic 
property is dynamically expressed by the determinantal weight, det [$? (Z k (T))], on paths 

l<j,k<N ? 

{Z k (t),te [0,T]}f =1 for any < T < oo. 

In the present paper, we would like to discuss the formula of Borodin and Corwin (1.10) 
for the O'Connell process from the view point of our theory of determinantal processes 
[16, 17, 14]. In order to do that, we first rewrite their expression. 



1.2 Main Results 

In the present paper, we set v = av = (az/^aP^, . . . ,aV^) e M. N and aV(-) = YljLi ^auj (•)• 
We first report the reexpression of the Fredholm determinant of Borodin and Corwin. 

Proposition 1.1 Assume that sup{|P ? | : 1 < j ' < N} < l/(2a) and {dj}^ =1 are all distinct. 
Fort> 0, (x,x') E R 2 , let 

N r 

K^ a (t;x,x') = J2 / dyp(t,x\V j )p(t,y\0)^ ,a ^ + ( L13 ) 

where 

F(l - a(r' - z)) V Y(al~-% (LU) 

r:P({r})=l, ^ ^ '' 



5 



and put 



K,n(x,x') = K.N(x,x';t,i/,a) 



Then, for fteR, 



= -Kf(l/t;x/t,x'/t) 



Det \S(x - x') - K, N (x,x')l {x/<h) 

(x,x')£R 2 L 
N ( -,\N' N ' rh 

Y~1^T] dx 3 det 



JC N (xj,x k ) 



(1.15) 



;i.i6) 



The points are following. 

(i) The expression of Borodin and Corwin can be rewritten as the Fredholm determinant 
of the rank iV operator with a kernel )Cn(x, x'), (x, x') G M 2 multiplied by an indicator 
l( x '<h)- Then the Fredholm series (1.16) has only N + 1 terms. 

(ii) The kernel K^ a (t; •, ■) is obtained from K^ N (t, -;t' , ■) of (1.1) by setting £(•) = z/(-) 
and t' = t and replacing the function $ p (-) by (•)• Equation (1.15) means that the 
kernel /Cat is the reciprocal-time transform of K^ a in the sense of [14]. 

Remark that an expression for the kernel /Cat, which is valid even when some of v/s 
coincide, is given by (3.5) in the proof of Proposition 1.1 in Section 3. It contains a contour 
integral on C(u) and, if {vj}f =1 are all distinct, the Cauchy integral is readily performed 
and (1.13) is obtained. For simplicity of expressions and arguments, here we assume that 
are all distinct. 

By the fact \im z ^ zT(z) = 1, in the combinatorial limit a — > 0, &p' a (z) — > §~(z), and 



thus 



\imK% a (t;x,x') = K%(t,x;t,x'), (x,x') G K 2 , t > 0, 



(1.17) 



where the rhs is the (equal time t! = t) correlation kernel (1.1) for the noncolliding BM 
without drift starting from a particle configuration z/(-). Then, the a — > limit of the rhs 
of (1.16) gives the Fredholm determinantal expression for the probability that all particle- 
positions are greater than the value ht, in the noncolliding BM without drift starting from 
v, when we observe the configuration at the reciprocal time 1/t; 



UmE au ' a 



Q a (x?(t) - h) = P^x^i/t) > ht], t>o,he 



(1.18) 



Note that our noncolliding Brownian motion, X{t) = (Xi(t), X 2 {t), . . . ,X N (t)), is ordered 
as Xi(t) < X 2 (t) < ■ ■ ■ < X^(t),t > in labeled configurations. By the reciprocal time 
relation proved in [14], the rhs of (1.18) is equal to the probability that all particle-positions 
are greater than h at time t, in the noncolliding BM with drift vector v, where all particles 
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are started from the origin. This initial state is given by the delta measure at the origin with 
multiplicity N, expressed by N5 , and we write the probability for this drifted noncolliding 
BM as Pp 5 °[ ■ }. Then we have 

(1.19) 



;i.2o) 



F'[Xi(l/i) > ht] = P£*°[Xi(t) >h], t>0,he 
Combining (1.18) and (1.19) gives the relation 



UmE au ' a 

a->0 



Q a (X?(t) - h) = P^*°[Xi(*) >h], t>0,he 



As mentioned before, 6 a (-) given by (1.9) is a geometric lifting of an indicator function 1(.> ). 
We will show that, in the combinatorial limit a — > 0, the transition probability density of 
the O'Connell process with av started according to the entrance law coming from "— oop" 
converges to that of the noncolliding BM with drift v started from N5q (see Lemma 2.1 in 
Section 2). Then, Proposition 1.1 will state that the result (1.10) by Borodin and Corwin 
[6] is a geometrical lifting of the Fredholm determinantal expression for the probability 
F^X^t) > h] of the drifted noncolliding BM. 

The complex function §~' a (z) appears in the kernel (1.13) for the O'Connell process is 
not entire; as shown by (1.14), it has simple poles at 



n 



z n = +r', neN= {1,2,3,...}. 

a 



1.21) 



(Note that all poles go to infinity in the limit a — > and the function becomes entire 
in the combinatorial limit.) Therefore, we will not obtain useful martingales to represent 
time evolutions of the system as in [17], but the single-time observables can have the CBM 
representations. The main result of the present paper is the following. 

For a configuration V(-) = J2f=i <%(') wrtn G M, 1 < j < A^, we consider the CBMs, 



Zj(t) starting from u j: 1 <j < N. That is, ^(0) 



and Wj(t) = 0, 1 < j < N. The 



expectation with respect to the CBMs under such an initial condition is denoted by 
Theorem 1.2 Under the same condition of Proposition 1.1, 







] = w 


det 




l<j,k<N . 



5 jk -<S>f a (Z k (l/t))l mm 



<ht) 



(1.22) 



The observable 6 a (Xf (t) — h), h G K is a softening of the indicator l(Xi(t)>/i). Theorem 
1.2 shows that its expectation for the O'Connell process started according to the entrance 
law coming from "— oop" has the determinantal CBM representation, in which the 'sharp' 
indicators l(y k (t)<ht)i 1 < k < N are observed, but the complex weight on paths is 'softened' 
as det [®$' a (Zk(-))], and the martingale property is lost. Further study of the maps 

of CBMs, $J' a (Z i (-)),l < j,k < N, and the system of CBMs with this determinantal 
weight is an interesting future problem. We hope that the present study will give some hint 
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for understanding why determinants appear in the processes which are not determinantal 
[6,7,1,8]. 

The paper is organized as follows. In Section 2 preliminaries of the O'Connell process 
and the noncolliding BM are given. The derivation of Proposition 1.1 from the result by 
Borodin and Corwin [6] is given in Section 3. Section 4 is devoted to the proof of Theorem 
1.2. Appendix A is prepared to give a sketch of a non-rigorous approach to deriving the 
result by Borodin and Corwin. (See the proof of Theorem 4.1.40 in [6] for a rigorous version 
of it.) 

2 O'Connell Process and Noncolliding Brownian Mo- 
tion 



2.1 Orthogonality and Recurrence Relations of Class-One Whit- 
taker Functions 

Let N — 2, 3, . . . and x = (xi, x 2 , ■ ■ ■ , x N ) G R N . The class-one Whittaker function ip^\x) 
has several integral representations, one of which was given by Givental [11], 



^\x)= [ exp(j^(T)W 

JTn(X) 



Here the integral is performed over the space Tn{x) of all real lower triangular arrays with 
size N,T = (T ijfc , 1 < k < j < N), with T N , k = x k , 1 < k < N, and 



N / j j-1 \ N-l j 



f$\t) = j>, J2 t ^ - E 3 ^.* -EE {e-^- T ^ + e-n + M-™y 



j=l \k=l k=l / j=l k=l 

We can prove that [22, 10] 



det [e x ^\ 

l^-^ix/a) = ^f N , (2.1) 
where h^iy) is the Vandermonde determinant 

M") = jt^n = n (^-^)- ( 2 - 2 ) 

~ h ~ l<j<t<N 

The following orthogonality relation is proved for the class-one Whittaker functions [25, 

29], 

I *!>)€'<^ = M^jvi S S(k ~ ° {k ' )] < (2 - 3) 
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for k, k' E R N , where sjv(-) is the density function of the Sklyanin measure [27] 



' l<j<£<N 



- ^M i nJfr--^> sinhT( r'" > [ <- R " ^ 

&n is the set of permutations of N indices, and er(fc') = {k'^y . . . , k' a ^) for a G &n- 
Borodin and Corwin proved that for a class of test functions, the orthogonality relation (2.3) 
can be extended for any k, k' e [6]. Moreover, the following recurrence relations with 
respect to v are established [18, 6]; for 1 < r < N — 1, v e C^, 

E II iT^^W^^f-E^Ww, < 2 - 5 > 

7C{1,...,AT}, jgj, V k 3 > \ j=l J 

\I\=r ke{l,2,...,N}\I 



where ej is the vector with ones in the slots of label I and zeros otherwise; 



1, 3 G /, 

0, je{l,...,7V}\7. 



In particular, for r = 1, 



JV 



j = l \<k<N:k^j V fe J ' 

where the £-th component of the vector e^y is (e{j})^ = 5^, 1 < j, £ < N. As fully discussed 
by Borodin and Corwin [6], the recurrence relations (2.5) are derived as the q — > 1 limit of 
the eigenfunction equations associated to the Macdonald difference operators in the theory 
of symmetric functions [19]. For more details on Whittaker functions, see [18, 2, 12, 22, 13, 6] 
and references therein. 



2.2 O'Connell Process 

In order to discuss the relationship between the O'Connell process and the noncolliding 
BM, we have introduced the parameter a > 0. The transition probability density for the 
O'Connell process with u is given by [14] 

P% a (t,y\x) = e -^lV^#M Q;(t)2/W) x,yeR N ,t> 0, (2.7) 

Wv (x/a) 

with 

Q a N (t,y\x) = j e^ 2 l^ N hx/a)i, {N) k {y/a)s N {ak)dk. (2.8) 
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(See also Proof of Proposition 4.1.32 in [6].) As a matter of fact, we can confirm that u(t, x) = 
P^ a (t,y\x) satisfies the Kolmogorov backward equation associated with the infinitesimal 
generator C u ^ a given by (1.6), 

l^ d 2 u(t,x) A d\og^\x/a) du(t,x) 

3=1 ■> J=l J J 

x E R N , t > 0, under the condition u(0, x) = 5{x — y) = n^LiO^.?' - 2/j)> 2/ e ^ JV - Assume that 
the initial configuration x e is given. Let Meff and < ii < £2 < • • • < *m < oo. Then, 
for this Markov process, the probability density function of the multi-time joint distributions 
is given by 



V^ a (t 1 , X M-t 2 ,xW-...-t M , X M) 

M-l 



m=l 



il>l } (x/a) m=1 

a-M G R^, 1 < m < M. 

Recall that we denote the O'Connell process by 

X a (t) = (X 1 a (t),X 2 a (t),...,X«(t)), t>0. (2.11) 

It is defined as an A^-particle diffusion process in K. such that its backward Kolmogorov 
equation is given by (2.9) and the finite-dimensional distributions are determined by (2.10). 
Therefore, (2.11) is a unique solution of the following stochastic differential equation for 
given initial configuration X a (0) = x e M. N , 



dX?{t) = dBj(t) + F^ a (X a {t)) dt, 1 < j < N,t > (2.12) 

with 



J 3 



F^ a (x) = V\og^\x/a), (2.13) 

where {Bj(t)}f =1 are independent one-dimensional standard BMs and [V]j denotes the j-th 
coordinate of a vector V. 

2.3 Special Entrance Law 

Let N eN, and define 

N-l N-l „ iV-1 , JV-1\ 
10 



O'Connell considered the process starting from x = —Mp and let M — > oo 
claimed in [22] (see also [2]) that 



-Mp) ~ Ce-^-^^exp (e M / 2 ^(T )) 



It was 



(2.15) 



as M — > oo, where the coefficient C and the critical point T° are independent of v. Then 
as a limit of (2.7) with (2.8), we have a probability density function 



P% a (t,x) = lim P£°(f,a:|-Mp) 

M— >oo 



with 



^(*,a!)= / e-*l fc l 2 /y^ fc (a;/a) Sw (afc)rffc 



(2.16) 



(2.17) 



for any t > 0. Since we have taken the limit M — >■ oo for the state —Mp, we cannot 
speak of initial configurations any longer, but for an arbitrary series of increasing times, 
< t 1 < t 2 < ■ ■ ■ < t M < oo, the probability density function of the multi-time joint 
distributions is given by 



M-l 



V% a (h, t 2 , x^- ...;t M , x™) = ] [ P^' a (Wi - ^, s^V" )?^, * (1) ) (2.18) 



m=l 



for e M. N , 1 < m < M. We can call the probability measure V^ a (t, x)dx with (2.16) and 



dx = Ylf=i dxj an entrance law coming from "— oop" [22] using a terminology of probability 
theory (see, for instance, Section XII. 4 of [24]). We note that, by (2.7), (2.18) is written as 

M-l 



rn=l 



The expectation with respect to the distribution of the present process started according 
to the special entrance law (2.16) is denoted by E"'°[-]. For measurable functions f( m \ 1 < 
m < M, 



M 



n/ M (* a (o) 

m=l 

e-^l^lW ITT f dx (m) I / (M) ( ^(M) ) ^) ( ^ ( M) /a)g ^ (tM _ tM _ i):r (M) |:r (M-l) ) 

U=V^ J 

M-l 

x ] [ /( m )(a ; ( m ))g^(t m -t m _ 1 , a; ( m )|x( m - 1 ))/( 1 )( : r( 1 ))^(t 1 , a; ( 1 )), (2.19) 



m=2 
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< h < ■ ■ ■ < t M < oo, where dx^ = ]jf =1 dx { ™\ 1 < m < M. 

Remark 1 The present special entrance law (2.16) is called a Whittaker measure by Borodin 
and Corwin [6] and denoted by WM^.^cc). Note that in the notation of [6], a Whittaker 
process is a 'triangular array extension' of the Whittaker measure and is not the same as 
the O'Connell process. 

When M = 1, for t > 0, (2.19) gives 
W , a[f{X a m 

= e -tmv^ f dxf{x) ^) {x / a)rN ^ x) 

Jr n 

= e -t\v\*/2a* f dxf{x)^\x/a) [ dke-^'^Ax/^s^ak). (2.20) 

Jrn J r n laK 

2.4 a — > Limit 

The transition probability density of the absorbing BM in W^r is given by the Karlin- 
McGregor determinant of (1.3), 

q N (t,y\x)= det \p(t, yj \x k )\, x,yeW N ,t>0. (2.21) 

l<j,k<N 

Consider the drift transform of (2.21), 

q v N (t,y\x) = exp |-^| 2 + v • (v - x) j q N (t, y\x). 

Then, if v G W N = {x G M N : x\ < x 2 < ■ ■ ■ < x N }, the transition probability density of 
the noncolliding BM with drift v is given by [3] 

det [e"' Vk ] 

p" N (t,y\x) = e ~ tm2/2 lJ ^J^M t, y\x), x,yeW N , t>0. (2.22) 

l<j,k<N 

In the limit Uj — > 0, 1 < j < N, of (2.22) the transition probability density of, the noncollid- 
ing BM is given by 

p N (t, y\x) = ^lq N (t, y\x), x,yeW N ,t> 0. (2.23) 
h N {x) 

We prove the following. (The superscript au is used for the processes with drift vector 
av.) 

Lemma 2.1 For v G Wat, 

XmiV a ^ a {t,x)dx = p N (r\x/t\is)d(x/t) 

= p" N (t,x\0)dx, t>0. (2.24) 
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Proof By the asymptotics (2.1) and the definition (2.21) of qw, we have 



l ima N(N-l)/2 e - m y 2 ^ a N) {x/c 



a-5>0 



t J 



h N (u) 



(2.25) 



For defined by the integral (2.17), we can show that the Whittaker function with 
purely imaginary index multiplied by the Sklyanin density, ijj^ N \(-)sN(ak), is uniformly 

— %dr\j 

integrable in a > with respect to the Gaussian measure e" -*l fe l 2 / 2 rffc,t > 0. Then the 
integral and the limit a — > is interchangeable. Since 



det [e~ ix ^ 



as a — > 



by (2.1), and (2.4) gives s N (ak) ~ a N< ^ N - 1 \h N (k)) 2 / {(2n) N N\} , as a -> 0, we have 
\ima- N( - N -^ 2 r N (t,x) 



, dke~ tlkl2/2 det \e- ix ^]h N (ik) 
(2n) N N\ J RN i<j/<N [ 1 NK ' 



+-N(N+l)/4 i 



(2^)^/2 " AT! 
By multi-linearity of determinant, 



d(Vik) det 

i<j,e<N 



-(v r tk i +ix j /Vi) 2 /2 



2/9 £-1 



2tt 



I J («Vi^ - iVtk m ) 



m=l 



1 

iV! 



d(Vtfc) det 

i<j,e<N 



-(Vtkt+ixj/Vt) 2 /^ 



2/0 £-1 



2tt 



(iVikg — i\ftk n 



= det 

i<j,e<N 



= det 

i<j,e<N 



d(Vtk)- 



m=l 



J 



2tt 



N2/o £-1 



Y[ {iVtk — iVtk v 



m=l 



e -(u+iXj/Vi) 2 /2 

du f= I I (iu — iVtk m ) 

m=l 



(2.26) 



The integral in the determinant (2.26) can be identified with an integral representation given 
by Bleher and Kuijlaars [5, 16] for the multiple Hermite polynomial of type II, 



p^ixj/Vt) with &-i(0 = EVtfc m (0- 



m=l 



(We set £o(") = an d rim=i(' '') = !•) It is a monic polynomial of Xj/ \ft with degree t — 1. 
Then (2.26) is equal to the Vandermonde determinant 

h N (x/Vt) =t N ( N ~ 1 V 4 h N (x/t). 
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Therefore, we obtain 



^a- N ^y^ N (t,x) = ^L^e-^h^x/t). (2.27) 

Combining (2.25) and (2.27), we obtain the equality 

lim^f,*) = ^j^- qN (t-\x/t\u)t- N , (2.28) 

which gives the first equality of (2.24) by the formula (2.23). The second equality is concluded 
by the reciprocal relation proved as Theorem 2.1 in [14]. The proof is then completed. | 



Remark 2 Moreover, if we take the limit v — > in (2.24), we have the following 

x\0) 

+-N 2 /2 

N — e-l*l 2 /2* (Ma;)) 2. (2.29) 



lim ]imV%' a (t,x) = p N (t,x\0) 



(2vr)^ 2 n-Lir(j) 

This is the probability density of the eigenvalue distribution of the Gaussian unitary ensemble 
(GUE) with variance a 2 = t of random matrix theory. It implies that a geometric lifting of 
the GUE-eigenvalue distribution is the v — >■ limit of the entrance law coming from "— oop" , 

V a N {t,x) = jjmptffrx) 

= 4 N \x/a)r N (t,x) 

= ^ N \x/a) [ e- t|fe|2/ V (A °fc(a;/a)sjv(afc)dfe- (2.30) 

J R N -™ K 



3 Proof of Proposition 1.1 

We start from the following result found as Theorem 4.1.40 in Borodin and Corwin [6]. 
(See Appendix A for a discussion of how this result relates to the properties of Whittaker 
functions as in (2.3) and (2.6).) Let 5 = sup{\uj\ : 1 < j ' < N} and choose < 5 < 1 so that 
5 < 5/2. Then for u G R 

oo L „ , 

E^[exp(-«e-^W)l = E 77 II f S , \K u (v v v k )}, (3.1) 

where K u (v, v') is given by (1.11). By assumption on 5, we can take the contour C(— is) such 
that any pair of v, v' G C(— v) satisfies \v — v'\ < 1. Then by using the identity 

1 

v + s — v' 



/ e- (v+s - v ' )b db, 
Jo 
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and interchanging the ds and db integrals, which is justified by giving appropriate decay 
bounds on the integrand, 



K u (v,v') = ^ dbe v ' b — r(-s)r(l + s)TT = T{V + V A , u > e -{v+')»*>'l*«*ft* 

JO J-i 



By multi-linearity of determinants, the rhs of (3.1) is equal to 

oo L 



L=0 j=l 



dv 3 > 



c(-u) Z-ni J 

ioo+5 



dbje v ' b ' 



x det 

l<j,k<L 



r{ Vj + u t ) 



J^oo+s^i 1 ; 1 nir( s + v j + v e ) 

^ ^ POD _ 

L=0 j=l 1/0 ~ 



oo ~ L 



with 



C(-l/) 



ioo+<5 



ds 



N 



T(v + u e ) 

l r(s + w + ^) 



Here we set 



XM s e -s6'+ti)s/a 2 +ts 2 /2a 2 -i)(6'-6) 



i/j = aPj, 1 < j < iV, it = e h/a , 
and change the integration variables in (3.2) and (3.3) to be 

bj — (h — Xj)/a, 1 < j ' < L, v = —aw, s = as. 

Then (3.2) is rewritten as 

I L rh 



L\ 

L=0 j=l 



with 



dw 



C(P) 
N 



2ni 



ioo+8 



d's 



K(x,x') = -a <p — — T(-as)T(l + as) 



ioo+8 



2m 



x r(a(^ - W)) c ( x '-tw)s+t^/2+w(x-x') 



A r(a(s + - it;)) 
where 5 = 5/a. Note that (3.5) is independent of /i. 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



15 



By assumption, {^}f=i are all distinct. Then the Cauchy integral with respect to w on 
Civ) is readily performed as follows. For each v j: l < j < N, 



Res 

w=Vj 

Since 



r(a(Pj- — w)) 
r(a(s + Vj — w)) J aT(as) 



r(-aS)r(l + as) 1 

= -1(1 — as) 



ar(as) a 
by residue calculation, (3.5) becomes 

N pioo+S 



K(x,x>) = -£/ ^r(l-as) 



x TT T(a(^ - Vj)) ( x '_tv i )8+fS 2 /2+v i (x-x') C3 g\ 

J- 1 r(a(s + P/-P,-)) ' 
i<e<N, v v 

Next, in each term of the summation over j, 1 < j < N, in (3.6), we change the integration 
variable, s~ — > y, as 

s = —(x'/t + iy) + Vj. 

Then (3.6) is written as 

K(x,x') = -J2 / . dj/r(l-a{P J --(x'/* + i?/)}) 

T(a(v e -Vj)) e x2 l 2t e -*(%-A) 2 /2 e -ty 2 /2 



n 



X ^.W-^/t + iy)})^ 



By definition of (1.3) 



-t(pj-x/t) 2 /2 ]_ g-V/2 1 

= _ ^r 1 ,^!^), — =- = -=p(r\y|o), t > o, 



and thus 

1 p x 2 /2t ^ /•oo+i(S+x7*-Pj) 

■ ' E ' 

X1 _^_ wt+W) n (3.7) 



]_ e x 2 /2t ^ r^+i(S+x' /t-Uj) 

K(x,x') = d yP {r\x/t\vj) P {r\ y \Q,) 
xI Yi uu , .-./m TT r («(^ ~ ^i)) 



1<^<JV: 



Here we consider each integral with respect to y in the summation. Note that p(t 1 ,y\0) 
and 1/T(a{v£ — (x'/t + iy)}), 1 < I < N, £ ^ j are all entire functions of y. The function 
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r(l — a{pj — (x r /t+iy)}) has simple poles, which are located at y n = i(n/a + x' /t — Uj),n G N. 
Since the assumption 5 < 1 gives 5 = 8 /a < 1/a, Qy n > 5 + x'/t — 9j,n G N, and thus the 
integrand has no singularity in the strip between the line C — {z — y+i(5+x' /t—Vj) :t/GR} 
and the real axis R in C, 1 < j < N. Owing to the Gaussian factor p(t~ x ,y\0), the integral 
on C can be replaced by that over R. Then we can conclude that 

K(x,x') = --JL^K**(l/t;x/t,zf/t), (3.8) 

where K^ a is given by (1.13). By the multi-linearity and the cyclic property (the gauge 
invariance) of determinants, det [K(xj, Xk)} = (— l) L t~ L det [K^ a (l/t; Xj/t, Xk/t)]. 

For fixed t > 0, a > 0, consider the integral operator in L 2 (R) with the kernel (1.13). It 
can be regarded as the projection on the subspace Span jp(t, : 1 < j ' < TV j, and has 

a domain given by Span j f R dyp(t, y\0)§p' a (- + iy) : 1 < j < iv|. As both subspaces have 
dimensions JV, det [Kj a (l/t, x,-/t, x fc /i)] = for L > N. Then (1.16) is valid and the 

l<j,k<L 

proof is completed. | 



4 Proof of Theorem 1.2 

Let x(') be a real integrable function and consider the following integral; for N' < N,t > 
0,a > 0, 

r N ' 

I NI [ X \ = / dx Hxix,) det m a (t;x 3 ,x k )}. (4.1) 

The determinant is defined using the notion of permutations and any permutation a E &n' 
can be decomposed into a product of cycles. Let the number of cycles in the decomposition 
be £(a) and express a by a = CiC 2 . . . Q( CT ). Here C\ denotes a cyclic permutation, 1 < A < 
£(cr), and if the size of a cycle is q\, it is written as C\ = (c\(l)c\(2) ■ ■ ■ c\(q\)), c\(j) G 
{1,2,..., N'}. By definition, we can assume the periodicity c x (j + q\) = c x (j), 1 < j < q x . 
Then 

t{tr) qx 

~ h ~ ves N , x=ij=i 

and (4.1) is written as 

Mxi= E(- i ) Ar ' _ ' (<T) n G ^x] 

o-GSjv/ A=l 

with 

G[c x , X ]= / I No?a; CA(j )x(a; CA ( i ))K^ a (t; 

) • (4-2) 
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Now we write (1.13) as 

K^ a (t;x,x') = [ D{dv) [ dyp(t,x\v)p(t,y\0)$l' a (x' + ty) (4.3) 

with P(-) = Ylf=i <%■(')> an< ^ rewrite (4.2) as 

C[c A ,x]= / ] J < ^c A (i)X(^ CA (i)) / ^(du CA (j)) 

x / ^c A (j+i)P(Mc A (j)k A (j^^ (4.4) 

Here note that, when we applied (4.3) to each 1 < j ' < q\, we labeled the integration variables 
as v — > f CA (j) and y — > y CA ( J+1 ) corresponding to x = x CA (j) an d a;' = x Cx (j +1 ), respectively. By 
Fubini's theorem, (4.4) is equal to 

\ (lHdv Cx(j) ) / \J\dx Cx{k) p(t,x Cx{k) \v Cx{k) )x(x Cx{ k))\ 
x / tl {^ A (^i)P(t,y CA (£+i)|0)$p AW ' a (a; CA ^ +1) + iy cx (e+i))\ 

n{x(K AW (i))^ (fc),a (^ A( , + i)W)} 



= w 

Then (4.1) becomes 



det 

l<j,k<N' 



K j ' a (Mt))xm)) 



By the Fredholm expansion formula for determinant, we obtain the equality 



N 



A"=0 



det 

l<j,k<N 



'6 jk -$f a (Z k (t)) X (V k (t)) 



(4.5) 



(4.6) 



By setting %(■) = 1(.<m, h G K. and performing the reciprocal time transform, the combina- 
tion of (1.16) and (4.6) gives (1.22). Then the proof is completed. ( 



Appendix 

A Sketch of a Non- Rigorous Derivation of (3.1) 

Here we provide a sketch of a non-rigorous approach to deriving Borodin and Corwin's 
Theorem 4.1.40 [6]. We work with the Whittaker functions and use their orthogonality (2.3) 
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and recurrence relations (2.6) to compute moments W' a [(e~ x *®/ a ) K ], k e N = {0, 1, 2, . . . }. 
This computation can be done rigorously as Section 4.1.4 of [6]. We then take a power 
series of these moments in an attempt to recover the Laplace transform of the distribution 
of e _x i^^ a , E^' a [exp(— ue~ x " (*)/°)], > 0. This is the place where the derivation becomes 
non-rigorous because the power series is divergent for all values of u and the moments do not 
identify the Laplace transform of the distribution. Nevertheless, proceeding formally and 
working with the divergent series we can recover the formula from Theorem 4.1.40 of [6]. 
Borodin and Corwin work at the higher level of g-Whittaker measures where the analogues 
of the moments are bounded by one and can be used to rigorously compute the g-deformed 
version of the Laplace transform of the distribution, which can be written as a Fredholm 
determinant. They then proved that the g-Whittaker measure converges weakly to the 
Whittaker measure, the g-deformed Laplace transform converges to the Laplace transform, 
and the Fredholm determinant has a limit which yields Theorem 4.1.40 of [6]. The fact that 
the formal calculations we describe actually recover the correct answer can be attributed to 
the fact that these are limits of the rigorous calculations done one level higher. 
The expectation at a single time t > given by (2.20) can be written as 

E»> a [f(X a (t))} 

-t\V\ 2 /2a 2 



= e 



dke- tlkl2/2 s N (ak) [ dxf(x)^\x/a)ij {N \(x/a). (A.l) 

J R N - mK 

First we consider the case with f(x) = e~ Xl / a . By (2.6), 

N 1 

e-*^(*/a) = £ n 77— 7<-Ue 0}) (*/«), (A.2) 



Ve — Vj 

j=l 1<£<N: t 3 



and 



/ dxe- Xl/a ^\x/a)^ N) h (x/a) 

J R N 



N _ 
j=l 1<£<N: * 3 JK 
N 1 1 

= a N y^ If , ttt7, ^ 5(ak - a(-iv + iei j \)), 

^ 11 v t - Vj s N (ak)N\ ^ m 

j = l 1<£<N: 1 3 ' ae&N 

where we used the orthogonality relation (2.3) extended to complex indices as in Section 
4.1.4 of [6]. Then (A.l) gives 

N 1 

E i/,o[ e -X 1 «(t)/o] = e -t\V\y2a 2 ^ Y[ ~ 



j=l 1<£<N: 
1*3 



vg. - Vj 



1 f t N 

x m S exp i (~ iu <p) + *( e mW))' 
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We can see 



1 

TV! 



<t£&N 

for v G R N . Then, if we set 



re6jv I p=i J 



/^(^) = e^ 2 n^ 7 , (A.3) 

£=1 1 



we have the expression 

E ,,a [e -X f (t)/a] = e t/2o» i —f/' a (v) 

Jc(-u) 2th 



C(-iz) 

t/2o»i" 1 /^>», f >o. (A.4) 



Next we consider (A.l) in the case f(x) = e 2x ^/ a . By (A. 2), 



1 



= E n ^f^-e^ 



j 1= l 1<£!<AT 



Applying the recurrence relation (2.6), it becomes 

N N J,( N ) 



Wle {n} -e {j JxM 



- - II II 



j 1= l l<j 2 <N: l<e!<N: ll n l<i 2 <N: V 2 n>y 31 



N 1 1 

+ E TT - J — TT — - — ^£W>/«)- 

ji=l l<^i<JV: " ■ ?1 l<fe<iV: 2 Jl 



Moreover, it is rewritten as 

N N u u 2 

E E — '1 - + l II II - ^(-iiz+ie^j+ie^jjC*/ ) 



l^v, — Z/o, + 1 - LJ - - LJ - Up — Uj 

j 1= l j 2 =l 02 11 r=1 i<^ r <7V: tr > 

N 



+e n n ue2+ \_ 4^e {J j*M- 

j 1= ll<4<A?": Jl £l 1<£ 2 <V * 2 - 71 
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Then, by using the orthogonality relation (2.3) and following the similar procedure to the 
first case, we have 

E ^,aj e -2X 1 a (t)/aj 



2t/2a 2 



N 

+ E 

J'i=l 



N N 



Ji=lj2=l 32 31 r=l 



e t(-u jr )/a 2 



n 

l<i r <N: 



t{-v n )/a 2 



n 



i<^i<w 



e *(-^'i+l)/a 2 



n 



— Z/f, + 1 + 



Here we consider a determinant of a matrix of size two 



det 

l<?'^<2 



1 




1 


f j + 1 — Vg 




1/(^2 + 1 -Ui) 












1 - (ui - w 2 ) 2 ' 



l/(ui + 1 - u 2 ) 



(A.5) 



which is equal to the symmetrization of (v 2 — i>i)/(i> 2 — fi+1) with respect to indices j G {1, 2} 

Of Vj% 



1)2 -V! V ± ~ V 2 



V 2 —Vi + \ V1—V2 + I 



Then we obtain the expression 



r 2 

i E ,,ar -2X 1 -( t )/ol = ^a 2 \ TT / ^ det 

2 L J 2^7 c( _, ) 27rii<^<2 



n 2 



1> 7 - + 1 — Vg 



r=l 



+ 



dv 



C(-i/) 2vrit; + 2 - w 



(A.6) 



By the similar calculation with the orthogonality relation (2.3) and the recurrence relation 
(2.6) of the Whittaker functions using the symmetrization identity 



e n 



Mi) ~ Mp) 



det 



+ 1 - ^ 



(A.7) 



one can prove the following. For any k£N 

_Lfl£">a|- e -«-X?(t)/aj — J_E^ a [( e - x i a W/ a ) K j 



^/2« 2 \p 1 TT 



A:|A| = 
i(A) 



det 



mi!m 2 ! • • • A -L J c( _^ 27U i<j,k<i(\) 



Vj + \j - w fe 



II {/^'>,)/^'>. + + \- - 1)}. (A.8) 

j'=i 
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where the summation is over all partitions 

A = (Ai, A 2 , . . . ) = l mi 2 m2 . . . , Ai > A 2 > • • • > 0, mj e N , j > 1 

conditioned that |A| = Ylj>i = K - Here /(A) denotes the length of A. (Precisely speaking, 
by using the orthogonality relation (2.3) and the recurrence relation (2.6) of the Whittaker 
functions, Borodin and Corwin gave a multiple contour-integral representation for general 
moment, W' a [{e~ x ^/ a )% k G N, in Lemma 4.1.29 in [6]. This integral formula involves 
nested contours. Then by deforming them to all be a small circle, denoted here by C(— v), 
Borodin and Corwin derived the formula (A. 8) (Proposition 6.2.7 in [6]), in which the identity 
(A. 7) was used.) 

One can rewrite (A. 8) in a suggestive form as 



±h e tiiM"^' 

L=0 n=(n u n 2 ,...,n L )m L : r=l J ^~ u ) 



x det 

l<j,k<L 



<i=l™i- 
1 



Vj + rij — Vk 



II + !)••• fx ■"(<■.! • n r 1)}. 



For 3?-u > one would like to recover the Laplace transform of the distribution of e x iW/ a 
from the moments via 



~^T^ a [(e-^^r] = E^ a [exp (- ue -^i W/a)] . 



K = 



One checks from (A. 8) that the moments grow super-exponentially, so this interchange of 
expectation and summation is unjustifiable and constitutes the physics 'replica trick'. Nev- 
ertheless we proceed now complete formally. By reordering terms in an unbounded manner, 
we arrive at the formula 



E u ' a [exp (-ue 



)] 



oo L 

Eif e n 

L=0 TleN L r=l 



x det 

i<i,fc<i 



C(-u) 
^rijt/20? 



dv r 
2tH 



i ")'■ f x J -"ir J )J7-"(r ! + !)••• ftf"( Vj + n 3 - 1) 



det 

2m i<j,k<L 



Vj + rij 

dV r 

'C(-u) 

which was given as (3.1) in the text, where 

e 



oo L 

E^n 

L=0 r=l 



K u {Vj,V k ) , 



(A.9) 



K u (v,v') = J2 



x nt/2a 2 



n=l 



v + n — v ■ 



7 (-u) n ftf*(v)f£**(v + !)••• fTiv + n - 1). (A.10) 
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By (A.3), 



r/'^r/'^v + 1) • • • r/> a (v + n -i) = n villi +\,y 

Then (A. 10) is equal to 

EJL r(u 4- uA ii n P tvn/a 2 +tn 2 /2a? 
(_l)nTT_L^±^^ . (A.ll) 

Since T(— s)T(l + s) = —71/ sin(7rs) by Euler's reflection formula and since — ir/ sin(7rs) has 
simple poles at s = n e Z with residues (— l) n , (A.ll) can be reexpressed as (1.11). 
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